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On a Random Matrix Models of Quantum Relaxation 

J. L. Lebowitz, A. Lytova, and L. Pastur 



Abstract. In paper [7] two of us (J.L. and L.P.) considered a matrix model 
for a two-level system interacting with a n X n reservoir and assuming that the 
interaction is modelled by a random matrix. We presented there a formula for 
the reduced density matrix in the limit n — > oo as well as several its properties 
and asymptotic forms in various regimes. In this paper we give the proofs of 
assertions, announced in \7\. We present also a new fact about the model (see 
Theorem I2.1| l as well as additional discussions of topics of [7] 



1. Introduction 

The model considered in [7] can be viewed as a random matrix version of 
the spin-boson model, widely used in studies of open quantum systems (see e.g. 
review works jS], I12j and references therein). We mention here that one of the first 
models of this type, namely the model where the classical system is represented by 
a harmonic oscillator coupled linearly with the oscillator reservoir, was considered 
by N. Bogolyubov in 1945 [3], Chapter IV. 

We recall now the model, proposed and discussed in [Tj. Let /i„ be a Hermitian 
nx n matrix with eigenvalues E^"\ j = 1, ...,n. We characterize the spectrum of 
hn by its normalized counting measure of eigenvalues 

(1.1) z.(")(A)=n-i5]xA(£;f^); / ^t\dE)^l, 

J=l 

where xa is the indicator of an interval A C R. We assume that z/q"'' converges 
weakly as n ^ c» to a limiting probability measure vq, i.e. that for any bounded 
and continuous function </? : M — > M we have: 

v[p\dE)'p{E) = / vo{dE)^{E), / v^idE) = 1. 

-OO J —OQ J —OO 

Let Wn be a Hermitian n x n random matrix, whose probability density is 
(1.3) Q-'cxp{-Ti-wl/2}, 
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where Qn is the normahzation constant. In other words, the entries Wjk, 1 < J < 
fc < n of the matrix Wn are independent Gaussian random variables with 

(1.4) {w,k) = 0, = 1, J, fc = 1, n, {{^w,kf) = {{^w,kf) 3 ¥^ K 

where the symbol (...) denotes here and below the expectation with respect to the 
distribution (jl.3p . This probability distribution is known as the Gaussian Unitary 
Ensemble (GUE) [9]. 

We define the Hamiltonian of our composite system S2,n as a random 2?! x 2n 
matrix of the form 

(1.5) = scr^ «) 1„ + I2 «) /i„ + va"" «) Wn/n^/^ 

where li (l = 2, n) is the / x / unit matrix, cr^ and are the Pauli matrices 

1 \ 2 f 1 

1 j''' =(,0 -1 

the symbol (gi denotes the tensor product, and s and v are positive parameters. 
The first term in ()1.5|) is the Hamiltonian of the two-level sytem ^2, the second 
term is the Hamiltonian of the 71- level system (reservoir) Sn , and the third term is 
an interaction between them. Thus s determines the energy scale of the isolated 
small system (2s is its level spacing), and v plays the role of the coupling constant 
between 52 and 5„. We write the Hamiltonian in the form 

(1.6) = ff^"' +Af("), 
where 

i?^") = sa' (g) 1„ + I2 «) hn, M^") = va'' (g) Wn/n^/^ 
and choose the basis in C^(8)C" in which the matrix Hq"^ is diagonal: 

(1.7) {H^J'^U^pk^X^^^^S^pS.k, \^:]=Ef+as, a,/3 = ±, j,k = \,..,n. 
Assume that at i = the density matrix of the composite system 52, n is 

(1.8) ^(:)(4"\o) = p(o)®Pfc, 

where p(0) is a 2 x 2 positive definite matrix of unit trace and Pk is the projection 
on the state of energy E^'^ of the reservoir. Let 11^"^ (t) be the density matrix of 
the composite system 52,n at time t, corresponding to the initial density matrix 
M^^O) of UHl): 

(1.9) M^")(i) = C/^"^(-t)/il"^(4"\0)f/(")(t), i7(")(i) = e^*^'"'. 
Then the reduced density matrix of the small system is defined as 

n 

i.e. p*^"^ is obtained from the density matrix (|1.9p of the whole composite (closed) 
system by tracing out the reservoir degrees of freedom. It follows from Theorem 
12.11 below that the variance of the reduced density matrix vanishes as n — > 00, i.e. 
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that selfaveraging. This aUows us to confine ourselves to the study of the 

mean reduced density matrix p*-"' (i?^"'' , t) : 

n 

(1.11) pS(4"',^)-E(^5a-(*)>= E T^a.(4"',^W7(o), 

J = l /3,7=± 

where 

n 

(1-12) = y: (uSp.i-t) + u^tsA*)) 

is the "transfer" matrix, an analog of the influence functional by Feynman- Vernon 

m- 

Notice that we can equally consider the factorized initial condition 0) in 

which the microcanonical distribution Pk of the reservoir is replaced by its canonical 
distribution e~^^" '' /Zg^\ We have evidently 

n 

(1.13) f,i-\p,0) = E^"'''^"VL"H4"\o)/4"^ 

It is also easy to write the corresponding reduced density matrix. 

2. Selfaveraging of Reduced Density Matrix 

Theorem 2.1. Let p^"'(i?^"\i) be the reduced density matrix il.lO\) of the 
composite system 5*2. „ — S2 + Sn, given by il.l]) - il.9\) . Then we have 

(2.1) Var|pi':l(4"\i)|<^, a,5 = ±. 

To prove the theorem we use the following facts. 

Proposition 2.2. (Poincare-Nash inequality). If a random Gaussian vector 
X — satisfies conditions {^j) — 0, (^j$^k) — Cjk, j, k — I, and functions 

$1,2 : ^ C have bounded partial derivatives, then 

(2.2) Cov{$i, $2} ((*!, $2)) - ($i> ($2) 

< ((CV$i,V$i))5((CV$2,V$2))^ 

where 

p 

(CV$,V$)= E ^jfe(V*)j(V*)fe- 

j,k=l 

For the proof of the inequality see e.g. [2], Theorem 1.6.4. 

Proposition 2.3. (Duhamel formula). If Mi, M2 are n x n-matrices, then 

(2.3) ^iM^+^h)t ^ gA/it ^ /■* eM^it-.)^j^^(M^ + M,).^^_ 

Jo 

The proof is elementary. Notice, that Duhamel formula allows us to obtain the 
derivative of the matrix {t) with respect to the entry wim, l,m = 1, n of the 
matrix Wn in (jl.Sp : 

(2 A) ^^Spkit) _n^r'y. („) _ („) 
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Proof, (of the Theorem \2.l\f . By using the Poincare-Nash inequaUty ()2.2p 
with 

n 

differentiation formula (|2.4p . and then Schwartz inequahty, we obtain 



(2.5) Var 5:C/S4.(0t/S,,(-0 

' TElf^i^Ki- ^)rE E^-i...(^)<U(-*) 



2i;2t 



n 



ds 



and here and below all the sums over the Latin indices will be from 1 to ti, and the 
sum over the Greek indices will be over ±. Notice that 



(2.6) 
and 



EI^-™./3.(-^)l =Elc^;t.(*-^)r-i' 



K,,l 



K.l 



(2-7) EEu':U-t+^)u\z,ii) =E(Et^5..K.H+.^)t^^"U(i-^) 

jl J2 K.l 



Hence, we have by (l23|), (EH), and ((2J)) : 



Var{Ef^S.(*)^SkH)}< 



4w2i2 



a, (5 = ±. 



Now, taking into account p.lOp and the fact that p(0) is a 2 x 2 positive definite 
matrix and of unit trace, we obtain (12.11). □ 



3. Equilibrium Properties 

We begin by considering the equilibrium (time independent) microcanonical 
density matrix of the composite system 52,„: 



(3.1) 



= 5{\ - i/(")) / Tr S{X - H(")) 



Following a standard prescription of statistical mechanics, we will replace the Dirac 
delta-function in (|3.ip by the function {2e)~^X£i where Xe is the indicator of the 
interval (—£,£), and e <C A. Then the reduced microcanonical density matrix, i.e. 
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the microcanonical density matrix of 52, n, traced with respect to the states of Sn, 
is the 2x2 matrix of the form 

(3.2) .(«)(A)- 



where 

n 

(3.3) T7(;)(A) = (2en)-i^Xs(A-iJ("))«,,,,. 

The corresponding canonical distribution of the composite system is 

(3.4) e-^""- /Tr e"^^" , 
and the reduced distribution of the small system is 

(3.5) ^—^ W : 

where (cf (1531) 1 

n 

(3.6) ^(")(A) = {z.(;)(A)}„,,=±, z.i!;)(A) = n-i^xA(i/("^)a,,^„ 

and XA is the indicator of an interval A of the spectral axis. 

Theorem 3.1. Consider the 2 x 2 matrix measure v^"^^ of hS.b]) . Then 
(i) there exists non-random diagonal 2x2 matrix measure 

such that the weak convergence: 

(3.7) lim = ly 



holds with probability 1; 

(a) if 



, A — Z 



(3.8) fM)=l T^-r^^^^^^^^^^ii-^>«-± 



is the Stieltjes transform ofva, and Vf) is defined by then the pair fa{z), a — 

± is a unique solution of the system of two coupled functional equations 

MdE) 
E + sa ~ z ~ v'^ f^a{z) ' 

in the class of functions analytic for 3z ^ 0, and satisfying the condition '^fa{z) ■ 
3z > 0, 3z ^ 0; 

(Hi) nonnegative measures Va, a = ± have the unit total mass, VaiM) — 1; 

and if the measure vq of hl.2\) is absolute continuous and sup < oo, then Ua, 

agr 

a = ± are also absolute continuous, and we have 

(3.9) u'J\)<s^x^l,'^{^i)■, 

(iv) for any A G M with probability 1 there exists the limit of the reduced micro- 
canonical distribution 

(3.10) lim = w, 
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where 

(3.11) c.(A)= ^^^^ 
and 

(3.12) I7q7 = 5a,7T'a, Va{X) = {2e)~'^ I'aidn), 

analogous formulas are also valid for the limits of the reduced canonical distribution 



Remark 3.2. The limiting measures Va, a = ± can be found from their 
Stiehjes transforms fa, a = ± via the inversion formula [Tl: 

(3.13) zyc«(A) = TT^Mim / Qfa{\ + iT)dX. 

To prove the theorem we need the following auxiliary fact. 

Proposition 3.3. Let $ be a function of n x n hermitian matrix, bounded 
together with its derivatives. Then we have for the GUE matrix Wn of (|1.3p : 

The proof of proposition follows from (|1.3|) - (|1.4|) and the integration by parts for- 
mula. 

Proof, (of the Theorem [3T|) . Denote 
(3.15) {z) = {H^''')~zy\^zj^Q 

the resolvent of (11.51) and set 



(3.16) gi"^{z) = n-'Y.GS,-,j(')- 

(n) 

It follows from the spectral theorem for Hermitian matrices that ga-y is the Stieltjes 
transform of and in view of the one-to-one correspondence between measures 
and their Stieltjes transforms (see [l], Section 59) to prove the weak convergence 
p.7p with probability 1 it suffices to prove that with probability 1 gi"' converges 
to Sajfa uniformly on a compact set of C\R. Denote 



(3.17) fL^H^) {9^^i^))-n^'T.iGi%i^))- 

For further purposes it is convenient to start by considering the functions 

n 

(3.18) <^(i) = «-'E(^S7.(^))' 

i=i 

where the matrix U^"\t) is defined in (|1.9p . By the spectral theorem for Hermitian 
matrices wi"'' is the Fourier transfrom of i^i"' and f^{z) is the generalized Fourier 
transform (see e.g. [llj) of u^aj(t): 

POO 

(3.19) f^a-}{z)=i-' e-'''ui"^{t)dt, -sz<0. 
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Notice that the matrix is diagonal with respect to the Latin indices. 

Indeed, since Wn in (II. 5|) is the GUE random matrix whose probabiUty law ()1.3|) 
is unitary invariant, we have for any unitary n x n-matrix U : 

(exp{iti/(")}) = {exp{it{H^"'^ + vn'^^'^G'' ® [/*)}). 

In particularly, for any diagonal unitary matrix U ~ {e"''^(5jfe}"j,^i with distinct 
i-Pj G [0, 27r), j = 1, n we obtain 

((exp{zti/(")})„,-^fe) = e*('^^-^'=)((exp{zti/(")})„,-^,). 

This implies 

(3-20) = C/^yO = (exp{*ti/(")}„,,^,). 

Hence we can write (|3.18D as 



1 . , 

(3.21) <H^) = -E<!.w- 

It follows now from jl^), (HJ]), and ((3ll)l that 



Hence taking into account (|3.20ll and (|3.2ip . we obtain: 

(3.23) V^S^ii)^ t dse^^-^^S 
Jo 

Jo Jo 



where 



rrf-EE(^-— iu=u-{u)). 

By using Schwartz inequality and inequality (|3.34p below we have that 
(3-24) \rilM<^- 

Here and below we use the notation C for all positive quantities that do not depend 
on n, t, z and indexes. 

We will use the notations g'""\z), /'■"-' (z) and R^^\z) for the generalized Fourier 
transforms (see p.l9p ) of the 2 x 2-matrices 

(3.25) C/j") = {t/S,.(0}o/3=±, u(-^ = W:^{t)U=±, rfit) = {r%mo.p=±- 
We have from the spectral theorem, ((3T^ . ((3:20)) . and ((3J5)) (cf (pl9j) ) 

(3.26) Gf\z)^r-' e'^^'uf{t)dt={{Gc.jm^))]^,p^^. 3z < 0. 
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This, ()3.17|) . and (|3.23|) lead to the matrix relation: 

+sa'-z- v^a-f(^\zy^)Gf\z) = 1 + Rf\z). 

Since the resolvent G*^"^ (z) possesses the property 

the matrix /^"^^z) possesses the same property G C^, and 

Thus the matrix E'j"-' + sa^ — z — v^a^ J^"^^ {z)<t^ is invertible, its inverse 

(3.27) f^'^\Ef\z) = (£;]"' + sa" - z - z;V"/(")(z)a")~^ 
admits the bound 

(3.28) ||/(")(i?]"\z)||<|3zr\ 
and equation for g')I^\z) takes the form 

(3.29) G;"^(z) = /(")(£;f\z) + /(")(i?f \z)i?5"^(z). 
Applying to the equation the operation we obtain in view of (|l.lp 

/oo 
i^(")(dii;)/(")(i?, z) + n-i ^ \ z)i?(")(z). 

Since the resolvent G*^"-' (z) is analytic if 7^ and bounded from above by |9z|~-'^, 
we have the bound 

(3.31) ||/(")(z)||<|5zr\ 

implying that the sequence {/^"•'(z)}„>i consists of functions, analytic and uni- 
formly bounded in n and in z hy rj^^ if |3z| > ryo > 0. Hence there exists analytic 
2x2 matrix function /(z), 3z ^ 0, such that ||/(z)|| < |9z|~^, and an infinite 
subsequence {/^"'°'' (^)}fc>i that converges to /(z) uniformly on any compact set 
of C \ R. This and estimates (I3.24p . p.28p allow us to pass to the limit Uk ^ 00 
in p.30p and obtain that the limit of any converging subsequence of the sequence 
{/^"^(z)}„>i satisfies the matrix functional equation 



(3.32) /(z)= / MdE)fiE,z), 
where 

(3.33) fiE,z) = iE + sa'~z^v'ay{z)a-)-\ \\f{E, z)\\ < \-sz\-\ 

The equation is uniquely soluble in the class of 2 x 2 matrix functions, analytic for 
3z 7^ 0, and such that 3z3(/(z)^,^) > 0, G C^. Indeed, for any two solutions 
/i, /2 of the class, and .9 = /i — /2 we have 

/oo 
MdE) [E + s<t' - z-wV^/i(z)cr^] "^vV^5(z)ct^ 
-00 

X [E + sa'- -z- wV^/2(z)ct^] ~\ 

and by (|1.2p . (|3.33p we obtain inequality ||(7(z)|| < w^|3z|~^||g(z)|| from which it 
follows that g{z) = for u5z < 1, hence for any 3z 7^ by analyticity. The 
solution of (|3.32p is diagonal, fa^ = fa^ap, and pair fa, a — zL satisfies system 
|. This follows from the unique solvability of (|3.8p . We can rewrite p.8p in 
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the form fa{z) = fa{z + f~a{z)), where /^(z) is the Stiehjes transform of the 
unit non- negative measure v^{E) ~ vq{E — as) . Since f^{z) possesses the property 
hm^^oo ?7|/a(«?y)l = 1 and \^{z + v'^ f-a{z))\ > then hm^^oo ^yl/ali??)! = 1 

and fa{z), a = ± are Stiehjes transforms of the unit non-negative measures Va{\) 
((3l3l) (see [l], Section 59). 

In addition, the Tchebyshev inequahty and bound (|3.35p below imply that for 
any e > 

n\&^{z) - g^-^{z)\ >e}< lvar{ff("^)(z)} < 
Hence the series 

oo 

5]P{|/(!;)(z)-g(';)(z)|>e} 

n=l 

converges for any e > 0, > ?7o > 0, and by the Borel-Cantelli lemma we have for 
any fixed z, > ??o > 0, lim„^oo yi""* (z) = /^^(z) with probability 1. With the 
same probability this limiting relation is valid for all points of an infinite countable 
sequence {zj}j>i, |9zj| > ryo > 0, possessing an accumulation point. Hence on 
any compact of C\R with probability 1 lim„^oo 9aj{z) = fa-y{z), and we have the 
weak convergence p.7p with the formulas p.l|) - (|3.8p . 

Let us prove assertion (ii) of theorem. It follows from the 

Q{z + v^U^{z))dE 



_^ {E + sa~ 3fJ(z + «2/-„(z)))2 + (3(z + «2/_„(z)))2 



5/q(z) < sup v'aip) 

= TTSUp Vq{^i). 

We have now by p.l3p 

i/„(A) < |A| sup 

This implies (ii). To prove (iii) we notice that by (ii) measures z/q, a = ± are 
continuous. Thus we can pass to the limit n cx) in (|3.3p . written as 



i")(A) = ^(")([A + £,A-e)). 



This and ([33) imply ((3T0)) - ((3TT|) . □ 

Lemma 3.4. Under the conditions of Theorem \3.1\ 



(3.34) Var{C/5^^,(t)} < Var{<)(t)} < 



(3.35) Var{G5)^,(z)} < Var{gH(z)} < 3z ^ 0. 

Proof. Acting as in the case of Theorem 12.11 we obtain (|3.34p . The differen- 
tiation formula for the resolvent 
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following from the resolvent identity, together with Poincare-Nash inequality ()2.2p 
imply the first inequality in (j3.35|) : 



n 

l.7n K. 



The second inequality in (|3.35p can be proved by a similar argument. □ 



We will return now to functions of variable t and find the n — > oo limits of the 
sequences {t/]"'(t)}„>i and - tT^ YJU C/j"^(0}n>i- 

Theorem 3.5. Consider the 2 -k 2 matrices {C^j"^ (i)}„>i and {u^"-\t)}n>i, 

defined in \3. 25\) and 113.21]) and choose a subsequence {E^^^} that converges to a 
given E of the support of vq of HI. 2]) . Then there exist the limits 

(3.36) Usit) = lim_ u{t) = Mm u^''\t), 



wh- 



ere 



i>N —irj 

(3.38) \\UEit)\\ = i yt>o, 



(3.37) Usit) ^ e"*f{E,z)dz:^^ \im / e''* f{E, z)dz, Vry > 0, 

^■^Jl 27r Af^oo 



and 

(3.39) u{t) = ^ J^dz e'^' J MdE)f{E,z) 
with f{E, z) defined m l[3J^) : 

(3.40) fc.p{E,z) = f^{E,z)5„p, faiE,z) = [E^-z-v\f-o.{z)T\ E„ = E+as. 

Proof. It follows from (|3.26p . (|3.29p . and the inversion formula for the genere- 
lized Fourier transform 11 that 

(3.41) C/j") [t) = g(") (£;]") ,t)+i-^ f Q(") (i?f) , t - s)rf {s)ds, 
where 

(3.42) Q(")(ij(;;),t) = -L f e'^'f^'-\E'fJ,z)dz^UE{t) 

+ i / ^) [E - E^;^ + i;V-(/(")(z) - f{z))a-] f{E, z)dz. 

The resolvent identity yields 



/(«)(£;("), z) = L. + ^ [sa^ + {E- Ef) + v^a-f(-\z)a-]f(-\Ef,zl 

and we have for sufficiently big r] = 

(3.43) ||/(")(£;(»),^)||<_A_, \\f{E,z)\\< 



z-E 
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This together with (I3.31|) allow us to pass to the limit under the integral in the 
r.h.s. of (|3.42|) and to show that it vanishes as rt ^ cxi. Moreover, we conclude 
that integral in the r.h.s. of (|3.42p is bounded uniformly in n and Vt > 0. The 
uniform boundedness of the matrix UEit) follows from the equalities {UE)apit) = 

{UE)aa{t)Sa/3 and 



Hence Q(")(£'j"\i) converges to UE{t) as 71 ^ oo and is uniformly bounded in n 

and t. This together with ([3^ . ((3?24ll and equality 1 1 C/j"'' (t) 1 1 = 1, Vt > give us 
(1^:77)1 and (13351). 

To prove (|3.39p notice first that we have from (j3.30p 

/> /•OC 

(3.44) = — / dze''* / v'^''\dE)f{E,z) 

27r Jl J-oo 

n nOO 

+ ^ Ldz I ^^"^ (dE) {E, z) - f{E, z)] 



+ ri / n-^y^Q(^\Ef\t-s)rf\s)ds. 
We integrate by parts with respect to z in the first integral to obtain in view of 

(loo)) 

- J- r .^:\dE) I dz e-^^±^^^4^. 

2^U-oo ° Jl {E^-z~v^f-4z)y 

It follows from (|332)) - ((333)) and dO]) that \\fa{z)\\ = \z\-'^{l + o{l)), \z\ oo and 
||/q(^)|| < Thus the integral with respect to z is bounded and continuous 

function of E. This and the weak convergence i/g"' to I'o (see ()1.2p ) yield the 
convergence of the first term of the r.h.s. of ()3.44|) to the r.h.s. of ()3.39|) . 
Furthermore, by using p.43p . (|3.3ip and p.ip we obtain 



oo 



\z\ 

<C\ [ u^^\dE)+ ( ^^ + , max ||/(")(y-^r;)-/(y- 777)11 

VJ\E\>T J\x\>AX^+T \y\<A+T 

For any e > choosing consequently A = A{e), T = T(e,A), No = No{e,A,T), 
and taking in account (|1.2p . and convergence f'^^^z) to f{z) on any compact set 
in C\ M, we obtain that the second term of the r.h.s. of p.44p vanishes as n — > oo. 
At last (|3?24)) yields for the third term of the r.h.s. of (I04l) : 

r^^— -Ell / '^-e-(*-)/(")(£;„z)|| 

-jy^^ r """"^ ^"^^^ l^\\UE{t~s)\\+ l^dz\\f^-\E,z)-f{E,z)\\j, 

and taking into account p.38p . (|3.43p . and (|3.3ip we conclude that the term also 
vanishes as 77 — * oo uniformly in t, varying on a compact interval. □ 
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4. Time Evolution 

We will prove now the main general result of [7], a formula for the limit as 
n ^ cxD of the expectation (|l.lip of the reduced density matrix (jl.iop of our model 
formula (4.7) of [t . 

Theorem 4.1. Consider the model of composite system, defined by lll.l\) - l[T73\) . 
Choose a subsequence of eigenvalues ofhn of il.5\) that converges to a certain 

E G supp vo ■ Then we have for the limit as n ~* oo of the expectation of the 

reduced density matrix Iil.l0\} uniformly in t varying on a finite interval: 

(4.1) p^AE,t):= lim (4'^) , i) ^ ^ / dz, [ dzie^*(---) 

^ fajE, Zi)fs{E, Z2)pa,s{0) + V^f-ajE, Zi)/-^(£', Z2) fg^sjzi, Z2) P-g-sjO) 
l-v'^faj{zi,Z2)f-a,-sizi,Z2) 

where Li — (—00 + ir^i, 00 + irji), L2 — (—00 — ii]2, 00 — M/2), rji > 0, 1J2 > 0; 

/OG 
MdE)fpiE,z,)f^{E,Z2) 
-00 

with fa{E, z) defined in {3.4-0^ . 

Proof. In view of (jl.lip it suffices to prove the following expression for the 
average transfer matrix (jl.l2p : 

(4.3) Tgp-,8{E, t) jim T^^f^^siEk^ ,t) 

where the "two-point" functions T^fi-ysiE, zi, Z2) are analytic in zi and in Z2 outside 
the real axis and have the form 

(4.4) fgpysiE,zi,Z2) = ff}{E,z^)f^{E,Z2) 

y-{5a,p5^^S + v'^S-a,/3S-y,5f-0,^^izi, Z2))[l - fp ^^{zi, Z2) f-0 -^{zi, ^2)]"^ 

Acting as in the proof of the Thcorem l3.1l i.e., by using the Duhamel formula (12. 3p . 
p.l4p and differentiation formula (12. 4p we obtain the relation (cf (|3.22p ) 



-v' rdse^^'^-^Kl f drY,T^:l_^^,{Eti^M.r)u%{s-r) 
Jo Jo ^ 

Jo Jo ^ 



+ / dse^^'^-^^''"lri%{t,,s), 

J 

where 

(4-6) ^S,(ii,i2) = -EE«-/3™Hi)t^SA(^2)), l4"^(ti,i2)|<l, 



n 

m— 1 J — 1 
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and 

(4.7) ri"^(ti,t2) 



t2 



+ r'^-EE(^-l„./..J--)^-?™...(*2)<_(r-ti)) 
It follows from Schwartz and Poincare-Nash inequalities and estimate p.34p that 



(4.8) 



Let T^^1^^(_E, zi, Z2), > 0, ^^2:2 < be generalized Fourier transform of T^^1,^(_E, ^1 



so that 



^i/375(^'*l'^2 



(27r) 



dz2 / d^ie'(*=^-*^^^)T^J .(£;,zi,Z2), 



where Li = (— cx) + 1771,00 + 4771), L2 = (—00 — 7772,00 — 7772), 771 > 0, 772 > 0. In 
view of relations 



we have 



I I diie-*^C/15,,(-ii) = GS,,(.i) 



a.j,[5k \ 



and (|4.5p yields 

(4.9) ?iy,.(4:\-i,-2) 
_ 1 



G^:l,sz.)5,s + E /i';l(-2)T(;i,,(4:), zi, Z2) 



E GL%,,„(.i)i^ri_,,,(.i, .2) + e,V(-i 



Z2) 



Here fT^") and 7^"^ are generalized Fourier transforms of X^"' and r^"^ of (|4.6|) and 
()4.7p respectively, and as it follows from ()4.6|) the absolute values of Kj^^^^g{zi, Z2) 
are bounded uniformly in 77 by |!3z2|~^. 

To write (|4.9p in the matrix form for any fixed pair a, (3 we denote K^^J , ^ 



(") 



r^'^'' the 2 x 2-matrices, which entries are {Kj^J)^s = K^J^g etc., an^ ^v^^ 
2 X 2-matrices with the entries {K~lf''')ys = K^^^ _^ ^, -y, (5 = ±, so 

^2) = /(")(4::\ Z2)?i5(^i, + ^2) 

X Gg,,Jzi)l2 +v'J2 G%^,Jzi)K-!,"Hzi,Z2) 
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Plugging expression for g["^(z2) we obtain 

(4.10) f(J(4:\zi,Z2) = i?S(<',zi,Z2) 

where reminder rI^^{E^^ , -^1,-^2) is a 2 x 2-matrix, and according to (I3.24p . (14.811 . 
and uniform boundedness of K^^\zi^ Z2) and f^''^\Em\ z), we have 

n 

(4.11) lim \0:^{Ei^^,z,,Z2)\\=O, lim ^ (£(,"), zi, z^)!! = 0. 

m— 1 

Applying the operation X!m^i (|4.10p with kn = m wc obtain: 

/oo 

/oo 
^(")(di?)/^;(i?, zi)/(")(i?, Z2) + n-^ Yl RSiE^:\zuZ2). 
-oo ™ 



This implies that for any fixed a, /3, 7, S the limiting values -^^0/375(^1,22) 
lim„^o ^'^p^si-^'^y -^2) and A'Q__/3^_7_i-(zi, Z2) satisfy the system of linear equatioi 

Kafj-fSizi, Z2) = fi3y{zi,Z2) [5ap5jS + v"^ Ka-p-j^Zi, Z2)] , 

where f0^{zi, Z2) are defined in (|4.2|) . Solving this system we obtain 

, . f-l3-^{zi,Z2)[5a^-p5-~f,5+v'^Sfj^~f{zi,Z2)5a,fjS^,s\ 

(4.1d) A„ _/3 _7,5(Zi,Z2) = ^ — -7 . 

1 - W^//3,7(^l' ^2)7-/3 -7(2^1,^2) 

Now we return to the variables ti, t2- It follows from (|4.10p and (|4.1ip that 

(4.14) T^p^siEMM)^ lim f dz2 f dz^ e'<-'^'^-''''^ 

X /^;'(4:\-i)/SH4!\-2) 

and we have to prove the equahty: 

(4.15) T^p^s{EM,t2)^7^ I dz2 f dz^e'^'-'^''^'^^ fp{E,z^)f^{E,Z2) 

X [SapS^S +v'^Ka-p.-^^sizi,Z2)], 



which together with (|4.13p yields (|4.3p . Notice that for any fixed non-real zi , Z2 the 
integrand of (|4.14p tends to integrand of (|4.15|) . but it has no an integrable majorant. 
Because of this fact we replace _^ ^(zi, Z2) in (|4.14p by the corresponding entry 
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of the r.h.s. matrix of (I4.12p to obtain 

X [/i";(M,^l)/i"yA^,^2)+«^ E /i"U(M,^l)£l,(A*,^2)i^a,«,.-.„.(zi,Z2)]. 

Here we denote J i>^^^\dii)dz2dzi — ^^^v'^\dii) Jj^^dz2 Jj^^dzi. Now it remains 
to prove that the foUowing expressions 

z.("^(dA.)dz2d^i zi) - f{E, zi)] Z2)/(")(M, ^l)/^"^^, Z2), 

v'^^'^ {dfi)dZ2dzif{E, Z,)f{E, Z2) Zl) - /(m, Z,)] /(") (m, Z2), 

v'^^\dll)dZ2dz^f{E, Zl)f{E, Z2)f{p, Z^)f{tl, Z2) [K^^\Z^,Z2) - K{ZI, Z2)] , 

K"^ ~ ^o] (c^Ai) / ci2^2 / dzif{E,zi)f{E,Z2)fifJ-,zi)f{ii,Z2)K{zi,Z2), 



where K^"\zi, Z2) — K^p-^si^^^ -^2), vanish as n ^ 00. 

Since ||/(")(^, zi)|| < r/f^ and there exist r/„ i = 1,2 such that Z2)|| < 

2|z2 - A*r\ \\f^"HEil\z^)\\ < 2\z,\-\ |3zj| > rj,, then the norm of the first 
expression is bounded by 

1/ 1 (.2||/H(,^)_/(,^)|| + |4^)_£;|) r.(")(d^)/ 

'/I Ai Nir 7lJz2||z2-/X| 

We have by Schwartz inequahty 



\dz2\ ^ f f°° dx f°° dx 



1/2 

TT 



lJz2 1 1^2-/^1 \J_^x'^+rj^ J_^{x- ^ly + rjjj 772' 

This and the uniform convergence of /'■"•' to / on a compact set of C \ K imply 
that the first expression vanishes as n — > 00. Treating similarly the remaining three 
expressions we prove that they tends to zero as n ^ cx3. □ 



5. Van-Hove Limit 

In this section we study the limiting case, where the coupling constant v of the 
system-reservoir interaction tends to zero, the time t tends to infinity while the 
transition rate, given by first order perturbation in the interaction, is kept fixed 
[H m [6l 110) . In terms of (14. 3|) this corresponds to making simultaneously the limits 

(5.1) V ^ 0, t ^ 00, T — tv^ fixed 

after the limit n — s- cx), i.e., in formula (|4.ip . 

We note that this limit as well as several other important topics of the small 
system-reservoir dynamics were considered by N.N. Bogolubov in 1945 [3] in the 
context of classical oscillator interacting linearly with the oscillator reservoir. 
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Theorem 5.1. Let the Fourier transform foiu) of the density Vq of the measure 
i>o in il.^) be absolutely integrable function: 



(5.2) / = Co < oo, 

J —oo 

/oo 
-oo 

Then the diagonal entries of the limiting reduced density matrix in in the van 

Hove limit are 



(5.4) pI"JE,t)=2'k 



where 

(5.5) T^{E)^2n[v'^{E) + v'^{E + 2as)]- 

and the off-diagonal entries are 

where fo is the Stiltjes transform of i/q : 

J— oo ^ ^ 

Lemma 5.2. In conditions 115. 115. 3\) of the Theorem \5.1\ next statements for 
the functions fa{z), a = ± are valid: 

(i) supQ2>o \fa{z)\ < Co, 

(ii) lim^^o i9/a(A + iO) = i/^(A - as), A G R. 

Proof. Estimate (i) follows from the representation of the functions fa{z), 
5z > in the form 

POO 

(5.7) /„(z) = i / e"'(-""+"+"'^-"("»^o(M)du 

Jo 

and condition 3z5/q,(z) > 0. It also follows from (|5.2p . (|5.7p that 



(5.8) lim/„(z) = /o(z-as) = W e'''^'-°"'> i7o{u)du, 3z > 0. 
'^^o Jo 

Hence 

(5.9) lim -3/a(A + iO) = -3? / e™(^""+^)i?o(u)dw = i^n(A - as), 
v^o n TT Jo 

Proof, (of the Theorem 15. ip . By using equalities (see (|4.2p ) 

/a. 0(^:1,^2) = Sz + V^Sf ' "^^ = ^1-^2, Sfa = faizi) ~ faiz2), 



□ 
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and by using analytic! ty of the integrand of (|4.ip in zi and in Z2, we can write 
the foUowing representation for the diagonal entries of the limiting reduced density 
matrix 



(5.10) 



PaAE,t)^^^^l dz2 I dz,e-'*''UE,z,)UE,Z2) 



LV, J LV 



^ Pa,a{0) f 1 ^^^^_^tg^V^Sfa5f-afa{E,Zi)fa{E,Z2) 



(27r)2 Sz{Sz + v^SU + v^Sf_o.) 

P-a,-«(0) f ^^f ^^^^_^,Sz^Hfc.{5z + vHf^)f.^{E,Z^)f.^{E,Z2) 



(27r)2 5z{5z + vHf^ + v^5f-^) 

^lUE,t)+lUE,t) + lUE,t), 

where LI — {zi : Szi = v^r/i}, L2 — {z2 : 5-22 = —v'^'112}, Vi ^^'^ V2 are arbitrarily 
chosen positive constants. 

To compute the limit (|5.ip of /{'(£', i) we change variables to — v^'^(zj — Ea), 
j — 1,2, and by Lemma l5.2l we have 



(2^)2 ''\,+f„(E + 2as + iO)JL, ^\2 + ME + 2as-iO) 

Computing last integrals by residues and applying equality 

(5.11) /o(A + iO) - /o(A - iO) = 27rzz/^(A) 
we obtain 

(5.12) vR-lim I^{E,t) = Pc.AO)e~^''^'''"^^+^"'\ 

where the symbol "vH-lim" denotes the double limit (|5.ip . 

Changing variables in I^{E, t) to C2 = v^'^{z2 ~ E^a) G ^2 = {C ^ Q'C = —^2}, 
Ci = v~'^{zi - Z2) e Li = {C : Q^C = ??i + m) yields 

(0) _..C, <5/a(Cl+'5/a) 



(5.13) /-(£;,0^ %r°i"^ / dC2 / dCie-^^^ 
(271-)^ Jl, 7l, 



(2^)2 -y^^ - Ci(Ci + '5,/c. + ^-a) 
1 1 



C2 + fo.{E-o. + v^C2) Ci + C2 + U{E^^ + v^iCi + C2)) ■ 

It follows from (15. 2p that the absolute value of integrand of (|5.14p is bounded from 
above by 



IC1IIC2IIC1+C2I ^Xi + (?7i + ry2)2\/AfT^V(Ai + A2)2 + r?^ 
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where c > does not depend on v, Xj — SftCji J — Now Schwartz inequahty 
yields for any B > 



Ai 



VAfTTv/(Ai + A2)2 + i 

dXi 



Ai 



dX. 



(A2 - + U/(A2 + + 1 



< 2 



dXo 



dXi 

Al Jo ((A2-^)2 + l)((A2 + ^)2 + l)i 



dAi 



dX. 



< oo. 



Al 7o ((A2 + ^)2 + l)i 
This aUows us to pass to hmit in integral (|5.13p by using (|5.8p and ( 15. lip 
vH- lim/^(£;, t) 

P-a.-a{0) 

Ai(Ai + 2TTi{v'Q{E) + v'^^{E - 2as)) 



— OO 

1 1 



A2 + fo{E - 2as - iO) Al + A2 + fo{E - 2as + iO)' 



Here integration path in Ai encircles zero from above. Computing last integrals by 
residues we have 



(5.14) vH-lim/3"(£:,i) = 2Trp^a~a{0) 



u'^{E-2as) ^_,r_^(E)''o{E-2as) 



T-c.{E) 



Treating similarly the term in the r.h.s. of (|5.10p we obtain 
(5.15) vH-lim/^(£;,t) 



27rpa,Q(0) 



rT^(E)^'o{E + 2as) 



Pa,a(0)e-^-^''o 



T-c.{E) 



TaiE) T^iE) 
Now the assertion (fO]) of theorem follows from the ([STTOj) . (l57T2|) . (ISH]) and ([STTS]) . 



dzie-''''U{E, zi)/_„(£;,02) 



Consider now the off-diagonal entry of (|4.1 



(5.16) Pc.,^c.{E,t) 



dzo 



dzy 



(2^)2 

it5^ v'^fa-a{zi,Z2)f-a,a{zi,Z2)fa{E, Zi)f^a{E, Z2) 



Pa,-a(0) 
(2^)2 



dzie 
dz2 



1 - l''*/a,-a(2:i, Z2)/-a,a(2:i, Z2) 



dzie 



-itSz • 



•'^fa,-aizi,Z2)f-a{E, Zi)fa{E, Z2) 
1 - V^fa-a{zi,Z2)f~a,a{zi., Z2) 



= '-^^^[ll{E,t)^Ili{E,t)]+Pj^^Il{E,t). 

To find the limit of Ii{E,t) of ()5.16p we change variables to (i ~ w^^(zi — Ea), 
C2 = v^'^{z2 - E^a)- This yields 



(5.17) 



I'iiE,t) = exp(-2asit) J^, 
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where 



(5.18) vH-limJl'(i?,i)= / fj^'^^' / 



and by 

II, C2 + fo{E ^2^7^70) Jl, Ci + fo{E + 2as + iQ) 
= (2Trf exp{ir((/o(i; + 2as + iO) ~ fo{E - 2as - iO))}. 

We have similarly: 

(5.19) = exp{-2asit) J^, 

„4 , 



TV _ f f „-irSC ^ fa-af-a,a 



1 1 



Cl+f-a.{Ea.+V^C,) C2+fa.{E-a.+V^C2)' 

Here we denote fp^^ — /^.^(-Eq + u^Cij-E'-a + ^"^(2) (see (|4.2I) '). By using the 
relations: 



«2(,5C + <5/-„,a)' ' ■ -4as + t;2(5C + <5/a,-«)' 
where (5//3^^ = //3(£^q + t'^Ci) - f-t{E-a + ^^(2), we obtain 

^' /— a.,a^ fa , — a 1 



(5.20) ^ f dC2 I dCie-'^"^-^ 

JL2 JLi ^ 



V'^fa,-af-a.a Cl ^ C2 + <5/-a,Q 
1 1 



Cl +/-«(£;« +1^2^l) C2+/a(£^-c«+t'2C2)' 

Notice that 



2co 
v'^{-qi + 772) 

Hence 



I /a -a I < 757:: Q; = ±. 



2co \' 1 



|1 - -f /a-a/a.-al > 1 " 1 \ 7 > , if ?7j > 2co, 

\{m+m)J 2 

and integrand of (|5.20p is uniformly bounded from above by integrable function 
C'dCi ~ C2||Ci|IC2|)^^- This allows us to pass to the limit in the integral in (|5.20p 
and obtain that 

(5.21) vH-lim J2''(£;,t) = 0. 
Treating similarly the term I^{E,t) of ()5.16|) we obtain 

(5.22) vH~lim/^'(^>0 = 0- 

Now assertion ()5.6p of the theorem following from ()5.16I1 - (I5.22|) . □ 



According to ()5.6p the off-diagonal entry of the reduced density matrix in the 
van Hove limit does not vanish but just oscillates as const-e"^*""*. The exponential 
that determines these fast oscillations (recall that t ^ oo) is the same as in the zero 
coupling (52-isolated) limit of our model (|1.5p . where the reduced density matrix is 
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hence is again const-e^^'"** ii a ^ /3, {a — —/3). 

In the case where the two-level system models a continuous quantum mechan- 
ical degree of freedom associated with a potential with two wells (see e.g. [8J for 
examples and discussion) , the above oscillation reflects the phase coherence between 
the quantum mechanical amplitudes for being in the left and right wells, a pure 
quantum mechanical effect. In this case our result means that an environment, 
modeled by a random matrix, does not destroy the quantum mechanical coherence, 
at least in the weak coupling regime corresponding to the van Hove limit. 

However, from the statistical mechanics point of view the absence of decay, 
moreover, fast oscillations, of the off-diagonal entries of the reduced density matrix 
seems not too natural. In this connection it worth noting that the fast ("micro- 
scopical") oscillating behaviour of pa/3, a ^ (3 can be converted into a decaying 
behaviour by several modification of our initial setting. 

One of them is to assume that the spacing 2s of our two-level system is random 
and continuously distributed, although concentrated around a certain 2so- In other 
words, it is necessary to assume that the two-level system is the subject of a certain 
(even small) noise. 

Another modification is to replace the van Hove limit 

(5.23) \im p{E,t)l2^r/t 

by 

pt+At 

(5.24) lim {2At)-^ p{E,ti) \^2^^u, dtt. 

t — ^oo,Ai — ^oo Jt — /\t 

If At = t, we just replace the limit t — > oo by the Cesaro limit (time average limit), 
a rather often used procedure in statistical mechanics. However the off-diagonal 
entry vanishes even for t oo, bat At/t 0, although with a smaller rate of 
decay. One can view this as an assumption on a sufficiently large (macroscopic) 
measurent time: « At « t. 
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